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Abstract: We construct a one parameter family of finite time blow ups to 
the co-rotational wave maps problem from S^ x R to S"^, parameterized by j^ € 
(i , f ] . The longitudinal function u{t, a) which is the main object of study will be 
obtained as a perturbation of a rescaled harmonic map of rotation index one from 
M^ to S*^. The domain of this harmonic map is identified with a neighborhood 
of the north pole in the domain 5^ via the exponential coordinates {a,0). In 
these coordinates u{t,a) — Q{X{t)a) +TZ{t,a), where Q{r) = 2arctanr, is the 
standard co-rotational harmonic map to the sphere, X{t) = t~^~'^, and Tl{t,a) 
is the error with local energy going to zero as i — )> 0. Blow up will occur at 
Vo'' (t, a) = (0, 0) due to energy concentration, and up to this point the solution will 

have regularity i?^"'''^". 
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1. Introduction 



2 , 1.1. Set-up and History. We study co-rotational wave maps [/ : S*^ x K — > 

5^ ^-> M"^ . Co- rotational means that the rotation index is one, or more precisely 

^^ ' U{t,ujx) = uiU{t,x) for {t,x) gRx S'^, and w € SO{2), acting in the standard 

way on S^. The metric on the target sphere is induced from the ambient M'^, and 

the domain is equipped with the standard metric of negative index 1, nip-y. With 

Ph ■ Einstein's summation convention in force, d^ = m^'^dj, irn^'') = {m/3-y)~^ , da 

■^ ' the associated volume element on the domain and (•, •) denoting the standard 

• , inner product on R^, wave maps are characterized by being critical with respect 

"•pi ' to the functional 

S ■ 

U^ I {dpU,d^U)da, ^ = 0,1,2. 

-Y-\ ' Recall that the energy 

>■ 
CN ; S{U)~ {DU{t,-),DU{t,-))dYo\s2 

is conserved for such critical points. If we use spherical coordinates a, 9 on the 
domain sphere and let u stand for the longitudinal angle on the target sphere, 
CO . we can think of U as {t,a,9) i— > (u(t, a),0), where u : (0,7r) — > (0, tt) is now a 

scalar map. u then satisfies 



sm(2M) 

— Utt + Uaa + cot au^ — 7^ — • (1 j 

2 sin a 

k> ' This equation will be the main object of study in this paper. The initial data 

; I , will be given at time tg small enough, and we will solve the equation backwards 

d • in time. Note that the problem at hand is critical in the energy sense. Shatah 

and Tahvildar-Zadeh [33] proved the existence of stationary solutions to this 
problem under less restrictive spatial equivariance assumptions. Their results 
include well-posedness of the Cauchy problem in the small energy setting. Here 
we are interested in constructing blow ups through an energy concentration 
scenario. Our construction will closely follow that of Krieger, Schlag and Tataru 



[Tl] where blow ups are constructed for co-rotational wave maps from K^ x R to 
S'^. The detailed discussion of the results in [M] is deferred to the next subsection. 
Here we briefly discuss some major recent developments, since the appearance 
of [13]. This is by no means an exhaustive account, although the reader can 
consult our references, in particular |14) and |26j . for a more comprehensive 
survey. Krieger et al. have used their method to construct blow ups for other 
wave type equations, in particular the NLW [IP and Yang-Mills [16] problems. 
In the setting of wave maps Carstea, |lj has extended the results of [H] to the 
case where the target is a more general surface of revolution. For the critical 
focusing wave equation, Donninger and Krieger [S] have recently modified these 
methods and constructed solutions which blow up at infinity. Going back to 
the wave map problem, Rodnianski and Sterbenz |27j considered fc— equivariant 
wave maps from IR^+^ to 5*^, for fc > 4. Using different techniques from those 
of [13], for each homotopy class k, they construct finite time blow up solutions 
which remain smooth until the blow up time. Moreover, they demonstrate that 
their blow up construction is stable under small equivariant perturbations of the 
initial data. In [26 Raphael and Rodnianski study the problem for all rotation 
numbers A; > 1 (as well as the 5*0(4) Yang-Mills), and are able to construct 
stable blow up solutions for all values of k. In their construction, for fc > 2 the 
prescribed stable blow up rate is independent of the choice of initial data in the 
open set of data leading to blow up, but there is an additional degree of freedom 
depending on the initial data for k = 1. More precisely, if T denotes the blow 
up time, for fc > 2 the blow up rate considered by Raphael and Rodnianski is 

A(t)=Cfc(l + o(l))- 



|log(T-t)|^ 
where as for fc > 1 they have 



A(i) = (T - i)e-Vli°g(^-*)l+0(i) . 

Related to this work, but in the context of Schrodinger equations, are the papers 
[i7],[18],[19],[20],[2T],[25],[22],[23]by Merle, Raphael, and Rodnianski, and 
a similar result by Perelman |24| . The construction of blow up wave maps by 
rescalings of harmonic maps are complemented by the result of Cote, Kenig, and 
Merle [3^ , where scattering is shown in the energy norm for the co-rotational case 
with initial data of energy smaller than the lowest energy harmonic map. See 
also [2] where the question of stability of harmonic map profiles is studied. We 
note that in the case where the dimension of the domain is larger than two, the 
study of singularities goes back to [4 and [5S] . As the focus here is the formation 
of singularities, we do not discuss the developments in the well-posedness and 
scattering directions in detail, but we refer the reader to the important papers 
[7]-[T3],[3n]-[3T], and [33]- [35], and the references therein. 



1.2. General Outline. After this short historical note we proceed to a general 
outline of the argument. As mentioned above, this outline will be very similar 
to that in 114] . It was shown by Struwe in |32] that if a co-rotational wave map 



from M^+^ to S'^ blows up, then a rescaled version of this map converges to a 
time-independent harmonic map from M^ to S'^ at the blow up point. Motivated 
by this, Krieger et al. imposed the ansatz u — Q + TZ where Q is a certain 
rescaling (to be specified below) of the harmonic map Q{r) = 2 arctan(r), and TZ 
is the remainder u — Q. They then find TZ in such a way that u becomes a wave 
map and the energy of TZ decays as i — ?► 0, leading to energy concentration at 
the origin. The idea here will be similar. The difference is that Q will no longer 
be a harmonic map from the whole space S*^ to the target 5^, but rather from 
R^ identified with a neighborhood of the north pole via spherical coordinates 
(a, 6). One can think of this M^ as the tangent space at the north pole, mapped 
to the sphere via the exponential map. More precisely we let 

u{t, a) = Q{X{t)a) + TZ{a, t), X{t) = i^^"". 
Here i' S (^, 1] is a fixed parameter. The reason for this choice of A is that ac- 



cording to [32] if we have blow up in the flat case, then there is a decomposition 
of u as above (with r denoting the radial coordinate instead of a) and a sequence 
of times ij — > with A(ii)|t,i| — ^ oo such that the rescaled functions u{ti, wj-t) 
converge to Q{r) in the strong energy topology. Note that even though this bub- 
bling result by Struwe - which to the best of my knowledge is not available in 
the case of a curved background - served as a motivation in |f4] , its validity was 
not used explicitly in the construction. 

The idea now is to rewrite ([T]) as an equation for TZ and use a fixed point 
argument relative to an appropriate norm to find TZ. Because we want the energy 
of TZ to decay as time goes to zero, our iteration norms will involve growth in 
time. However, the error coming from the failure of the rescaled harmonic map 
Q{X{t)a) to be an exact wave map, is too large to be placed in these iteration 
spaces (at least the ones used in [14] and this paper). For this reason we have 
to further breakdown the error TZ as Q{X{t)a) + TZ{ta) — U2k-i{t, a) + e(i, a), 
where U2k-i is an approximate solution to the wave maps equation, and e is 
the error playing the role that TZ played before. The index 2k — 1 indicates 
that the approximate solution will be refined iteratively, starting with UQ{t, a) — 
Q{X{t)a). The point is that if k is sufficiently large, the error due to M2fc-i not 
being an exact wave map is small enough to be placed in our iteration spaces. 
Of course we will need to make sure that the energy of U2k-i — Q decays as we 
approach the blow up time. Before we state the main theorem we remind the 
reader that in the equivariant setting the energy is given, up to a constant, by 



£{u) 



'^ 'l/.,2 , 2^ , sin^iu) 



[Uf +u„ 



2 sin^ a 



sin ada. 



and the local energy relative to the origin is defined as 

sin ada. 



^ / -^ /* [1/ 2 2^ sin^(u) 
-2 2sm a 



Theorem 1. Let v £ (21 1] ^^ arbitrary and to > sufficiently small. Define 
X{t) ~ t^^^" and fix a large integer N. Then there is a function u'^ satisfying\}\ 



^ Following 1141 we call u'^ the elliptic modifier. Also for noninteger 13, C^ means C^Pi'P IpI. 



u" eC''+^^-{{0<t<to, r<t}), 

fiocK) W < (t(A(t)))-2| logi|2, ast^O, 

and a blow up solution u to (Op in (0, to] which is of the form 

u{t, a) = Q{X{t)a) + u''{t, a) + e{t, a), < a < t, 
where e decays as t — > 0; more precisely, 

set^H^+^-iS^), etet^-^H''-{S^), £i,,ie){t)<t^, as i ^ 0, 

with spatial norms which are uniformly controlled as t ^>- 0. The solution 
u{t,a) extends as an H^^^^ solution to all of (0,to] x S"^, and its energy con- 
centrates in the cuspidal region < r < -w^y leading to blow up at r = t = 0. 



1.3. Overview of the Proof. In this subsection we provide a more detailed ac- 
count of the blow up construction and the proof of theorem [T] The proof here 
will combine the method of proof of the corresponding theorem in |14j with a 
perturbative argument. By this we mean that equation ([l} is replaced by its flat 
analogue 

1 sin(2M) 

-Utt + Uaa + —Ua ^ n , 

and the difference is included in the source term. In the appendix, we have 
collected the results from [T3] which can be used directly in our setting. 



Section 2 is devoted to constructing the approximate solution which we refer 
to as the elliptic modifier. The construction consists of adding successive refine- 
ments to the previous approximate solutions, starting with uo(t, a) = Q{t,a). 
This is done in two steps: first near a = and then near the boundary of the 
cone a — t. The error at step 2fc — 1 

/ «2 , o2 , . Q ^ Sin(2u2fc-l) 

2 sm a 
will satisfy 

'i?(log(2 + i?,))2'=-i- 



62*:-! = O 



t-2fe+2i2(iA)2 



Here R — Xa is the rescaled radial variable, and the 0{-) terms are uniform 
in r < t. < t < to. Note that for i/ > 1 the decay in time is slower than the 
decay of the corresponding error e2k-i in the flat case from (llj which satisfies 



e2fe-i=o( 



fl(log(2 + R)f''-^ 
i2(tA)2fc 



Even though we restrict to the case v < 1, the renormahzation step is carried 
out for all values oi v > ^ . This slower decay is the result of the error coming 
from replacing the curved wave equation by its flat counterpart, and is observed 
in the analysis near the boundary of the cone. It entails that if we could take 
v > 1, to be able to close the final fixed point argument we would need more it- 
erations of our approximation scheme (see section 5). In our analysis we will use 
the variable R near a = and the self-similar variable a = j near the boundary 
of the cone a — 1. The approximate solution U2k-i and the error e2k-i will 
involve expressions of the form (1 — a)"" 2 log™(l — a), which dictate the upper 
bound on the regularity of our final solution. In both cases we will encounter a 
Sturm-Liouville ODE which is why we use the name "elliptic modifier." Even 
though the arguments in this section are similar to those in [2] we are not able 
to use most of the results from that paper directly. In particular to be able to 
track the error from our perturbative argument we will need to introduce extra 
auxiliary variables in our function spaces (see section 2 for precise definitions), 
and we will use a slightly watered down version of the function algebras Q and 
Q' (in particular we do not track the exact number of log(l — a) factors at each 
iteration, and hence we dispense with the index on the algebras Qk and QJ^,). 



Section 3 corresponds to sections 4 and 7 in [Tl] and does not contain new 
results. We make the ansatz u{t,a) = U2k-i{t,a) +e{t,a) and recast equation 
P^ as an equation for e as 

1 cos(2Q(\a)) ,^ , , -rr / X 

-£tt+£aa + -£a n £ = N2k-l(£) + Hi[e) + e2k-l- 

a a^ 

Here N2k-i contains the nonlinearity coming from the original equation on 
the curved background and is defined in ^, and Hi, which is responsible for 
the perturbative error coming from replacing the curved equation with the flat 
one, is deflned in ©. We will change variables to i? = Aa and r = ^t~'^ (so 
now the blow up time corresponds to t = cx)) and rewrite the equation above in 
terms of the new variable e(T, R) = R'^e{t{T),a{T, R)) as 



'^^ T-t'^ \ -1 / ^T \ 2 -1 o / '^^ ~ 



= A-2i?^(iV2fc-l(i?-^£)+i?l(i?-5£)+e2fe-l), 



where 

r . .2 , 3 8 

"^ AR^ (l + i?2)2- 

£ is a self adjoint operator on L^((0, 00), dR) (see the appendix for the defi- 
nition of the domain of L) and its spectral properties, which are analyzed in [5] 
and [2] are summarized in the appendix. In particular with p{^)d^ denoting the 
spectral measure associated with C we consider the distorted Fourier transform 

■^ : / ^ /(C) = lim / 4>{R,0fiR)dR, 



with inverse 

:F-^ ■■ f^fir) - lim [%{R,0f{0p{0dt 
We apply T to the equation satisfied by e and let x = J-e to get 



-{dr-2^d^yx^cx 



A7- j^ / r) ^ Ar ^ r^ \ A_ / 2 ^ 



2^/c(a. - 2^e9c)^ + ^(/C^ - Y^/C + 2[/C,^ad): 



Here the operator /C is defined as 

Rdii.u = —2£_d^u + ICu, 

and as such encodes the deviation of ROr from being diagonalizable in the 
Fourier basis. We denote by H the fundamental solution operator of the operator 
on the left hand side of the equation for x. The mapping properties of JC and H 
are summarized in the appendix. 



In section 4 we introduce the Sobolev norms associated with the operator 
C. The final fixed point argument is carried out on the spaces defined by these 
norms. On the frequency side we define the weighted L^ norms 



/iL?^^:-( ri/(ono'v(eR)' 



/o 
For functions of the spacial variable R we define 

To control the decay in time we introduce the /,^^°o2j^'* spaces with norm 

L^°°Hp is defined similarly. The norm of the N2k-i term arising in the 
equation satisfied by the Fourier coefficient x is controlled in [14 , and here is 
recorded in the appendix. The new ingredient in the nonlinearity is ifi, and it 
is this nonlinearity which is responsible for the upper bound 1/ < 1. In particu- 
lar, in the renormalization step (Section 2) we carry out our computations for 
general v > \- 



Finally in section 5 we assemble the results from the previous sections to find 
the Fourier coefficient x = Te using a contraction mapping argument. We will 
use certain L°°'^Lp'* norms for this fixed point argument, and we need to take 
fc large enough depending on N to be able to place the term involving 62^-1 in 
our iteration space. We will then retrace our steps back to find e and verify that 
it satisfies the appropriate energy decay as i — ^ 0. 
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2. The Elliptic Profile Modifier 

We will construct approximate solutions to the wave maps equation in the light 
cone Co at the origin 



Co := {{t, a) : 0<a<t, <t< t^}. 

The approximations will be perturbations of a time-dependent harmonic map 
profile 

uo{t,a)=Q{X{t)a), \{t) ^ r^-". 

We will closely follow the construction in [14 , while making the necessary 
modifications necessary to handle the curved metric on the background. To be- 
gin, we introduce the following set of auxiliary variables: 

'^-(l°^(2 + ^^»'=t-(log(2 + i?^)f. 



u\ - 


{txr 


hr, - 


(log(2 + i?2))2 


f^2 ~ 


(rA)2 


h- 


= «', 


64 = 


= (tA)-2 ^ t^- 


65 = 


^{e\)-^ = t' 


R- 


- Aa, 




a 


a = 


' 1' 



Technically speaking, we do not need 63 because 63 = 0^65. However, as even 
analytic functions of a arise naturally in our computations, we prefer to keep this 
as a separate variable. Since we are working in the light cone Cq, the variables bi 
will be restricted to small intervals [0,6o,j]i respectively, and a will take values 
in [0,1]. Let 



r? := [0, 1] X [0, (») X [0, 6o,i] X • • • X [0, 60.5] • 

and denote the projection of J? onto the last six factors by Qa^ and the projec- 
tion onto the last five factors by f2a.R- We will work with asymptotic expansions 
of the functions we are concerned with in order to study their growth and decay. 
To be able to make optimal use of the new set of variables we need a few defini- 
tions. These definitions, and the statement of theorem[5]are all that will be used 
from this section in the remainder of the paper. The reader may therefore wish 
to skip the proof of the theorem in the first reading. The definitions below are 
technical. For the sake of readability the reader may think oi I S"^ [R'^ (log RY , Q) 
as the space of functions satisfying the following properties. They are analytic 
near R = and vanish to order m there. For large R they behave roughly like 
i?'^(logi?)™. Q means that near the boundary of the cone a = 1, the functions 
may contain singularities of the form (1 — a)''+2(log(l — a))" for some integer 
n. If Q is replaced by Q' , the exponent t^ + ■^ should be replaced by j^ — ^- The 
auxiliary variables bi are harmless and may be safely ignored. We now provide 
the precise definitions. 



Definition 1. Q is the algebra of continuous functions q : [0, 1] — > M with the 
following properties: 

(i) q is analytic in [0, 1) with an even expansion at 0. 

(ii) Near a = 1 we have an absolutely convergent expansion of the form, 

i — 00 CX3 

q = qo{a) + ^ ((^ " «)^''-'^''+^ E 92.-i..(«)(log(l - a)y 

i=l j=0 

00 

with analytic coefficients qo, qij, such that for each i only finitely many of 
the qi_j are not identically equal to zero. 

Definition 2. Q' is the space of continuous functions q : [0, 1] ^ M with the 
following properties: 

(i) q is analytic in [0, 1) with an even expansion at 0. 

(ii) Near a = 1 we have an absolutely convergent expansion of the form 

i — 00 OQ 

q = qo{a) + ^ ((^ " ^^'''^"'^ E 'Z2.-i.,(a)(log(l ~ a)y 

i=l j=0 

00 

(l-a)2-+E92..,(a)(log(l-a)F") 

with analytic coefficients qa, qi_j, such that for each i only finitely many of 
the qij are not identically equal to zero. 



Definition 3. S'™(i?'^(logi?)') is the class of analytic functions v : [0,oo) -^ M. 
with the following properties: 

(i) V vanishes of order m at R = and v(R) = i?™ '^''i=o^ ^j^^'' /"'^ small R. 
(a) V has a convergent expansion near R — go, 

v^ J2 c,,R''-^\\ogRy 

0<j<l+i 

Finally, 

Definition 4. a) S'™(i?'^(log Kf ,0) is the class of analytic functions f : J7 — J> M 
so that 

(i) V is analytic as a function of R, hi 

V : f2a-^ Q 
(ii) V vanishes to order m at R ^ and is of the form 

j = ao 

i;«i?™^Cj(a,&i,62,&3,&4,&5)i?'^' 

3=0 

around R ~ 0. 

(Hi) V has a convergent expansion at R — oo, 

v{R,-M,h2MMM)= J2 c,,{-MMMMM)R'''^\\ogRy 

0<j<l+l 

where the coefficients Cij : Qa.R -^ Q o,fe analytic with respect to hi. 

b) I S"'^{R^ {log Ry,Q) is the class of analytic functions w on the cone Cq 
which can be represented as 

w{t, a) = v{R, a, h^M, &3, &4, ^s), v € S-^ {R^ {\og R)\ Q). 

Note that this last representation is in general not unique. Note also that 
these definitions are slightly different from the corresponding ones in [14] . 

Our goal in this section is to prove the following theorem. 
Theorem 2. For any fc G N, there exists an approximate solution 

U2k-i e Q{R) + -^/5^(i?logi?, Q) 



to the wave maps equation ([ip, so that the corresponding error e2k-i satisfi, 



es 



fe-i 



It is worth mentioning that for i/ > 1, the decay in time of the error here 
|-^2fc+2i/-4'j jg glower than that of the corresponding error in [T3] (i^*^^"^). 



10 

Proof. Our strategy will be the same as that in [TJ- We start with uq as a 
first approximation and improve our approximations inductively by adding a 
the correction vt at the fc'th step: 

Uk = Uk~l +Vk- 

The error at this step is 

/ „, „, \ sin(2ufc) 

-Gfe = - c/j + o a + cot aOa W/t 5 — • 

V / 2 sin a 

If u were an exact solution, then the difference e = u — Uk-i would satisfy 

C0s(2Mfc_i) . 

^9^ TT-^ sm(2e) 

2 sm a 

sm{2uk-i) , , .. , , 

+ ^ . 2 (l-cos(2e)) = efe_i, (2) 

2sm a 

where Dg = —d^+d^a + cotada- As in [T3] we linearize this equation around 
e = 0, and substitute uq for Uk-i to get 

cos(2mo)\ 

Dg- o • 2 Jg-efc-i- (3) 

2 sm a ' 

For r ^ t the time derivative should play a less important role and we 
approximate the equation above by 



02 , .0 cos(2uo)\ 
Ojj + cot aOa r^ I £ = 



= efe-i, 
sm" a 

In fact we make one further simplification and replace sin a by a and cos a 
byl: 

ff^2 , lo cos(2'Uo)\ 

It is understood that the error ek-i includes the one previously defined and 
the new error resulting from this last simplification (this will be written out 
completely below). This simplification allows us to use the results from [T3] 
more directly. Similarly for r w t, where we will later use the self similar variable 
a, we approximate cos(2uo) by 1 and set sin a « a, cos a ~ 1, and rewrite ^ 
as 

(-df + dl + -do^-^)e = ek-i. 
\ a a^ / 

The heuristic arguments above lead us to the following two step improvement 
of our approximations. First we consider the region r <^ t and define V2k+i as a 
solution of 

„2 1 ^ cos(2mo)\ n />x 

with zero initial data at a = (to be interpreted suitably below, as the 
coefficients are singular at r = 0). Here e^^, is the "large" part of the error e2k 



11 



and will be defined during the iteration. In the second step we improve our 
approximation by adding V2k+2 defined by 



(^-d^ + dl + ^do.-^) V2k+2 = e§,+i . (5) 

Again the zero initial data at r = needs proper interpretation, and Cjj. i j^ is 
the principal part of the error e2k+i, to be defined below. Comparing with ^ 
we compute the successive errors as 



e2k = e2fe-i + N2k{v2k) + Ho{v2k), 

e2k+i = e\k - dt'V2k+i + N2k+i{v2k+i) + Hi{v2k^ 

Here 



cos(2u2fe) - cos(2uo) sin(2u2fc), .„ x ,x 

N2k+i{v) = r^ V + . „ — (cos(2i;) - 1) 

sm a 2 sm a 

+ -7r—2 (sin(2f ) - 2w) , (6) 

2 sm a 



respectively 



cos(2u2fe-i) - 1 , sin(2u2fc-i). ,^ -. ,x 

N2k{v) = r^ V -\ -— ^ (cos(2u) - 1) 

sm a z sm a 

+ '-^f^isH2v)-2v), (7) 

2sm a 



and 



HQ{v)=(--cota)do,v+(—-—^-)v, (8) 



respectively 



„. . /I , \„ , /cos 2wo cos(2mo)\ 

^i(f^) = cota 9qU+ 5 r^ )v. (9) 



12 



To carry out the construction just outlined and prove the theorem, we imple- 
ment the following four step induction scheme. For each fc > 1, 



fc-i 



V2k-i e ^ 



1 



^ (t''A)2j(tA)2(fe 



-IS'{R{\ogRf''-\Q) 



k-l 



^'^-- e E 77^7TW77IW^^^^(^(l°g^)""^ 2') 



i=o 



(i'^A)2j(tA)2(fe-j) 



fc-l 



^- ^ E 77^7TW77OTTI^^^'(^'(l°g^)""^ 2) 



fc-i 



(i''A)2j(iA)2(fe+i-j) 



t'ezfc e ^ 



^L. (i.A)2..(iA)2('c-.) 



[I S\R-' {log Rf^Q) 



J2b,IS\R{\ogRr-\Q') 



(10) 

(11) 
(12) 



(13) 



StepO: t^eo £ /S'^i?"^)- 



q2 a2 1 Q , sin(2uo) 

+ ( cot a)daUa + (— ^-^ —^ 

a 2 sm a 2a^ 



) sin(2uo) 



52(g(Aa)) + (- - cota)9„(Q(Aa)) + (— ^^ - -i^) sin(2Q(Aa)). 
a 2 sm a 2a^ 



So, 



t'eo 



(1 + ^^) 



4i? 



1^(1 + iy)- 



2R 



+ t' 



(l + i?2)2 ^ a+i?2 

2(sin2a-a2)\ ( R{1 - R^) 



a^ 



^sin'a J\il + R^)^ 



sin a — a cos a \ R 



o? sin a / 1 + i?2 



Note that the functions of a in the second term are analytic and even, and 

5- 



that i^ = 65. It follows that 



t^eo e /S'i(i?-i), 

as desired. 

Step 1 : define V2h-\ such that ITUl) holds. 

Define the principal part e2;._2 of e2fe-2 by setting 6^ = 0, i = 1 • • • 5, in (J13p 
with fc replaced by /c — 1. We can write e\y._2 = J2j=o ^2k~2 7' with 
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The non principal part t^e\i^_2 ~ t'^i^2k-2 — ^2k-2) belongs to 

^T. ^,.,)2,'ark-, IS\RilogRr-\Q'), 

which can be included in e2k-i- The inclusion above is clear for i 7^ 3, and 
follows by writing 63 ~ a^ /{t'^Xf' for i = 3. 

Now as in (|4]), the idea is to define V2k-i.j by requiring that 
/^2 , In cos(2uo) 



Note that if v[R) is a function of only R then {d^ + ^9^ - £25g^)w = A^Lw, 
where 

2,1. cos(2uo) _ 1 l"6fl^ + i?^ 

L_a« + -5H-— ^^^-a^ + -9«- ^2(1 + ^2)2- (14) 

Using this as a motivation we define V2k-i.j by requiring that 

[t\fLv2k-i,,=t^elk-2.y 

In doing so, we view the equation above as an ODE in R and treat all other 
variables as parameters. The errors resulting from this simplification will be 
studied in the next step. According to lemma [1] from the appendix V2k-i,j G 
if^xyHlx)2i>.-.^ IS^iRi^ogR)^'^-\ Q). Therefore if we define 

fc-2 

V2k-1 := /^V2k-lJ, 
3=0 

we have V2k-i e E,to ^t^x)^,^ly,.-n ISHR{logR)^'^-\ Q)0 

Step 2 : show that with V2k-i defined as in the previous step, (jlip holds. 

Thinking of V2k~i as a function of t, R, and a we write 

e2k-i = elk-2 + N2k~i{v2k-i) + Hi{v2k-i) + E*V2k-i + E°'V2k-i- 

2 In fact we have V2k-i S Ej=o (f^A)2j(t\)^(''-i) ^'^'^^^^'"S^^^'''^' ^^^ '^''* "^"^ '^° """^ '''''' 
this. 
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e\j._2 ^^^ taken care of in the previous step. N2k-i{v2k-i) and Hi{v2k-i) 
account for the contributions from the nonhnearity and approximating sin a by 
a. They are given by ([5]) and Q respectively. E^V2k-i contains the terms in 
dt'^2k-i{t, R,o-) where no derivatives apply to a and E°'V2k~i the terms in 

(^~df+d^,+^dr)v2k-l{t,R,a) 

where at least one derivative applies to a. The only new contribution com- 
pared with [2] is Hi{v2k-i), which we now bound. 

/cos(2uo) cos(2uo)\ sin^ a — a^ 



\ sm a — a 

]V2k-l = 5"^-^ COS(2uo)W2fc 

/ a'' sm a 



a^ sin a ^ a-^ sm' a 

fe-i 



■^g („,p.(L)-.>-.. ^"'W'°«"'""'g')- 



because the coefRcient is an even analytic function of a and cos(2wo) € 'S'''(l)- 
For the other term in Ii\ we write 

----E(^i;^^PIf^' ^.^^s\mo^^r-\Q\ (15) 

and compute 



f COta \daV2k-\ 



, a cos a — sin a s ^^ adaWj (a, R) + ROrWj (a, R) 



a) "" ^"'^ ^ a2sina ^^ (i''A)2j(tA)2(fe-J) 

Since the coefficient is an even analytic function of a and since 

this last expression can be placed in J2jZo — it^xviitxyif'-'i) • "^^^^ takes 
care of Hi{v2k-i) ■ The contributions of the other three terms can be analyzed 
just as in [Tl], but we reproduce the calculations for completeness, starting with 
N2k-i{v2k~i) ■ First note that adding up the Vj we get 

U2k-2-u„^j^IS\R\ogR,Q). (16) 

For this note that extra powers of l/(iA)2 (respectively l/^f^Xf) can be 
written as powers of 64 (respectively 65), and as such can safely be placed in 
IS^{1). Also note that from lemma [2] in the appendix, if 

v^j^IS\R\ogR,Q), 
then 
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sin«G -^IS^{R\ogR,Q), cosv £ IS%1,Q). 

{txy 

Using this and (fTO)) we compute 

cos(2mo) — cos(2u2fc-2) = 2cos(2uo)sin {u2k-2 ~ "o) 

+ 2 sin(2Mo) sin(M2fe-2 - "o) cos(u2fc-2 - uq) 
1 



JS\R\\ogR)\Q) 
1 



(U) 

IS\\ogR,Q) 



(tX) 
Hence, 



2COS(2-Uo) — COs(2u2fc-2) Ct^ {tX)^{cOs{2uo) — COs(2u2A;-2)) 
t r"2 'V2k-l — ^ f^ V2k-1 

sm a sm a R 



X 



E77^7TW7^YW^^^'(^(l«g^)""^2) 



j=o 



(t''A)2j(tA)2(fe-J 



^1 (/55(i?-l(l0gi^)2^Q) + ^/57(^(iog^)2fe+l^ g) 

X 7" rTTr^T~TTTm ^^ 



-^ (i''A)2j(U)2('=-j) 

fe-i 



(i'^A)2j(a)2(fe-j) 

^E (,.,)2.(a)2(.-, ^^H^aog^)"-^go 



i=o 



Here in the step before the last we have pulled out a factor of b^ from the 
second factor, and given away a factor of R^ to gain a factor of log R in the first 
factor. For the next term in N2k-i{v2k-i), note also that 

1 - COs(2w2fc-l) = vlk-i9{v2k-2) 

for some analytic function g. Moreover writing l/(tA)2 = 64 and l/(i''A)2 = 
65, we can think of V2k-i as belonging to /S''^(i?(logi?)2'^^i, Q) whenever we do 
not need the extra time decay. Therefore since g has a convergent Taylor series 
at 0, g(v2k-2) e 5'"(1, Q), and hence 

k — l 

i-cos(2.2.-i) e (E (,.;^)2,-(,\)2(.-,) ^^'(^(l»g^)'''^Q 
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Using this observation we compute 

2sin(2u2fe_2) .-, ,„ ^^ 

t ^ . 2 (1 - cos(2i;2fe-i)) 

zsm a 






(i'^A)2i(tA)2('^-0 



4=0 
fc-1 



^.^jt''A)2.(tA)2(fc-.r 

We have again pulled out factors oi b\b™, m + I = k, to pass to the last 
inclusion. Finally, using a similar observation as above, 

2COs(2u2/c-2) .„ . .„ .. 

t r^ (2t;2fe-i -sm(2w2fc-i)) 

sm a 



-IS' 



'a^2)(E 



- ^2 V ,-^;^^ (i'^A)2j(tA)2(fe-j) 

fc-1 

^ E fi.A)2a+/+™)(iA)2(3fc-,-i-™)-2^^'Wl0g^)''"''2) 
j,l,m=0 ^ ' ^ ' 

k — 1 

^E ft.,)2.-(a)2(.-.) ^^^wiogi^r-\s). 

This concludes the analysis of A'^2fc-i(''^2fe-i)- For E^V2k~i we can ignore the 
dependence on a and therefore the observation that 



^fe-l rri^/T^/n r^N9I._1\ fc-1 

2 



/^3(i?(l0gfl)2fc-l) x ^^ /^l(i?(l0gfi)2fc~l) 

^ (i'^A)2j(iA)2('=-j) J - ^ (i"^A)2j(a)2('=-j) ' 



shows that the contribution of this term can be placed in the right space. 
Finally for E°'V2k~i we use the representation (jlSp again, together with the fact 
that 

ada, a-^ da, il- 0^)91 -Q^Q'- (17) 

The contribution of (^tu^r(f'x^2(k-n to t^E''v2k-i is 
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+ (1 - a')dlW, - 2adaWj\ 
In view of (J17p , if we sum the expression above over j we see that 

as desired. 

Step 3 : define V2k such that (fT2)) holds. 

Write e2fc-i = 2^;^q 62^-1,/ with t^e2k-i,i & — (f-A)^'(«>)^'''"'' ' ^*^**^^g 
6i = we write the asymptotic component f2k-i,i of the principal part of e2fe-i,z 
near R = 00 as 



R ^''^^ 



(t''A)2'(iA)2('=-0 



3=0 

2fc-l 



T Yl aqjia){\ogRy, 



with qj G Q'. In view of ([S]) we consider 



t^{-df +dl + -da - ^)w2k.l = t^f2k-l, 

and homogeneity considerations lead us to seek solutions of the form 



2fc-l 

(t''A)2'(U)2('=-0- 



^2fe,/ = u.^^9J.u^^2(k-l)-l Yl ^2fe.;(a)(log i?F 



j=0 



Matching the powers of log R we see that W|^ ^ has to satisfy 



<'(-3f + aJ + i<j.-J,)(- 



W^ 



2k,i \ _ agj(g) - Fj{a) 



a a-/ v(t''A)2'(U)2('=-0-iy (^^A)2'(tA)2('=-0-i ^ 

where with /3 = /J(fc, I) =21 + v{2{k -I)- 1), 
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F,{a)^{j + l)[{{l + ^)l3 + a-')W, 



i+i 

2k,l 



{j + 2)ij + l)ia-^-{l + ,.r)Wi+f. (18) 

We are using the convention T4^|j, ^ = for j > 2fc here, and will solve this 
system successively in j as j decreases from 2fc — 1 to 0. Conjugating the power 
of t we get 

t'(~ {dt + ^-f + 52 + U^ - ^Mlu.i = «9.(«) - FM- 

With 

Lfi -^ (1 - a^)dl + {a-^ + 2a/3 - 2a)9a + {-f + /3 - a^^)^ (^g) 

we rewrite this last equation as 

L0Wi^,^aq,{a)~F,[a). (20) 

Our goal is to show that with zero Cauchy data imposed at a = 0, 

W^lfc,ie«^2' j = 0,l,-- -,2/0-1, ? = 0,l,...,fc-l. (21) 

For this we will use lemma |3] in the appendix. In particular, we use the fact 
that for /3 > 1/2 the unique solution to the equation 

Lpw = /, w;(0) = 0, daw{^) = 0, 
has the form 

w{a) = Ci0i(a) + C2<j)2{a) + C30i(a) / (f>2{a')qi{a')f{a')da' 



+ C402(a) / 0i(a')9i(a')/(a')rfa'- (22) 



Here oq < 1 is some number close to 1 , the Ci are constants which may depend 
on /3, and the fundamental solutions (pi can be written as 

Ma)=ri{a), Ma) = r2ia)il - a)^+K (23) 

with ri, r2 analytic, if /? — ^ ^ Z+, and (j)i has to be modified to 

(f>i{a)^ri{a)+c(j)2{a)\og{l~a), (24) 

for some constant c (depending on /3) if /3— i G Z+ . The function qi appearing 
in the particular solution can be written as 

gi(a) = (l-a)-^-^r3(a), (25) 
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with r^ analytic near a = 1. 

The following integral identity, valid for 7 7^ —1 and positive integers m, will 
also be useful 



1 

{l-a'y[\og{l-a')Y''da' 



^ -[log(l-a)r(l-ar+i 



7 + 1 
- ^ /'(I - «')"[log(l - a'r-'da'. (26) 

We will prove (PT|) inductively on k. We will write c for generic constants, 
which may depend on (3 and v. Note that near a = 1, proving (|2ip for fixed I is 
tantamount to finding an absolutely convergent expansion 



W^,/a) = qo{a) + ^^((l - a)(2-i'-^+^ ^ fe-i.,(a)(log(l - a))^" 

1=1 j=0 

+ (1 - a)2-+i ^ 92., (a)(log(l - a))^-) , (27) 

such that 50 and qij satisfy the conditions of definition [TJ To start off the 
induction, note that when fc = 1, / = so we are in the same setting as in |14] . 
and (PU)) becomes 

Li^W2^Q{a) = ah{a), 

L.Wloia) = (c - a~^)Wlo{a) + {ca - a~')daWlo{a), 

where h is analytic with an even expansion at a = 0. The behavior near a = 
is a direct consequence of part (i) of lemma [3] For a close to 1 we denote by hi 
and gi generic functions that are analytic near a = 1, and use (|22p - (|26p to write 

Wl,{a) = go{a) + g,{a){l - ar+i , 

W%{a) - ho{a) + h,(a){l - 0)-^+^ + h2{a){l - aY+i log(l - a), 

a ly ~ ^ ^ Z+ , and 

Wloia) = go(a) +5i(a)(l - a)''+H 52(a)(l - a)''+5 log(l - a), 

2 
Wl„ - h,{a) + (1 - ar+^ Y. h,n+i{a)[\og{l - a)]" 

m=0 
2 

+ (l-a)2''+i^;i™+4(a)[log(l-a)r, 



m—0 
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\i v~ ^ e Z+ . This finishes the first step of the induction. Now we assume we 
have proved (j2ip up to fc — 1 and prove it for k. Here we are assuming that once 
we estabUsh (PT|) we can prove ([T^ and p^ as well, and therefore we assume 
that qj e Q' in (PO]) . That this can be done, is shown below in the remainder of 
step 3 and in step 4. 

We fix / < fc — 1 , and write /3 for 

p{k,l)=^2l + {2{k-l)-l)v. (28) 

We will generically write gj, hj, gij or hij for functions that are analytic 
near a ~ 1. The first equation we need to solve is 



^/3^2M ^aq2k-i{a). 

Again the behavior near a = is a consequence of part (i) of lemma |3l For 
the behavior bear a = 1, note that the solution is given by the right hand side 
of (|22p with /(a) = aq2k-i{a)- The contribution of Ci0i + C2(t>2 can be written 
as 

go{a) + gi(a)(l - af+-- + 52(a)(l - a)''+^ log(l - a). 
Note that if / = then 

and otherwise 

where 34(1) — 0, and m is an integer. This shows that the contribution of 
Ci0i + C202 has the right form. Noting that (j)2qi is analytic near a = 1, we can 
write the contribution of (pi J (j)2qifda' as 



h go{a) + E((l - «)^'^-'^'^+^ E .92.-i,,(«)(log(l - a)y 

00 \ 

+ (l-a)2-+i^.g2,,,(a)(log(l-a))^ 



Regardless of whether j/ is rational, this can be written as 



''' This is not circular; in fact one can view this step as part of the larger iteration scheme 
lfT0ll-|[T3t. 
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OO 

OO 

+ (1 - a)(2-i)-+^+i ^ /,2,^, (a)(log(l - a))^" 

OO 

+ (1 - a)2-+^+i E /i2.-i.,(a)(log(l - a)y) . 

3=0 

According to (l28)) . 

for some integer m > i (depending on i), and 

(1 _ a)2-+^+i = (1 - a)(2"-i)''+5 (1 - a) 

for some integer n > i + 1 (depending on i), and therefore the expression 
above can be included in the right hand side of (|27l) (that for fixed i the sums 
over J are finite follows from the corresponding fact for /). We consider the 
contribution of (l>2 J (jjiqifda' next. Here we may get terms of the form 

c(l-a)-i[log(l-a)]™ 

in the integrand, which cannot be dealt with using (P^ . and which contribute 
extra logarithms. To make this more precise, write 

(i)i{a) = .9o(a) + gi(a)(l - af^^ log(l - a). 

First we consider the contribution of the analytic part g^. This gives integrals 
of the forms 



h^{a){l ~ a)"+'^ / g2,-i,,(a')(l-«')^''"'''"''"'(log(l-«')rrf«', 

J ao 

/io(a)(l - a)^+5 r g2,,,(a')(l - a')'^''-^-^(log(l - a')Yda'. 



ao 



Both of these have the right form, and only add an extra power of log(l — a) 
if (2z — 1)1/ — /3 — 1 or 2iiy — /? — ^ are negative integers. For the contribution of 
gi{a){l — a)'^+2 log(l — a), note that 

(?i(a).gi(a)(l - af+i log(l - a) ^ hi{a) log(l - a), 
so we need to consider integrals of the form 



22 



ho{a){l - af+i / ff2.-ij(a')(l - a')<'^''-'^^-^log{l - a')y+'da\ or 



which again have the right form. This completes the case oi j — 2k — 1. 

For j < 2k ~ 1, we can assume by induction that W^2fc i ^^'^ ^2ki have the 
right form, and in particular they have a cubic Taylor expansion at a = 0. It 
follows from psp that aqj (a) — Fj (a) has a Taylor expansion at a = beginning 
with a linear term. The behavior of M^jfe ; at a = is therefore a consequence of 
part (i) of lemma [3l Moreover, Q C Q' implies that aqj{a) — Fj{a) S Q', and 
the expansion of VFjj. ; follows in exactly the same way as above. This finishes 
the proof of (^1]) . 

We cannot use W2k,i for V2k,u because logi? is singular at i? = 0. However, in 
view of the computations above, we define 



2fe-l 

i 

V2kJ 



2k-l _. 



(rA)2'(a)2(fe-o .^^ 

and analyze the error near R = generated by this modification in the next 
step. Pulhng out a factor of a^ = R^/{tX)^ we see that 



Step 4 : show that with V2k defined as in the previous step, ([T^ holds. 
We write 



t'^e2k =i^(e2fe-i - e^k-i) 

a a^ 

+ t^N2kiv2k)+t^Hoiv2k), 

with N2k and Hq as in ([7]) and ^ respectively, and 
The first term has the form 
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fc-l 



i'(e2fc-i-e°;,_i)e^ — 



^^ (rA)2'(iA)2('=-0 



2fe /-i'\ 



/5^ (i?-^ (log i^)"^Q 



^M^i(i?(logi?)2'=^\Q'] 



The last sum is contained (J13p . For the first term we can write, with w € 
IS\R-\\ogR)^'',Q'), 



w = {1 — a^)w + 






e IS\R-\\ogRf\ Q) + 61/51 (i?(logi?)2'=-i, Q') 



and therefore the first sum can be placed in ([T^ too. According to the com- 
putations in step 3, the second term in e2k would be zero if we had log R instead 
of ^ log(l + i?2) in both 62^, and V2k- The difference is therefore obtained, if we 
replace the derivatives of i log(l + R"^) by derivatives of log R in 



a 



2k — 1 

Ew^4,(«)(^iog(i+i?2) 



3=0 

This difference can be written as a sum of expressions of the form 



(i'^A)2'(tA)2(fc-0 



2fe-l 

tE 

j=o 



wLM) 



5°(i?-2)(log(l+i?2)y-l 
+ 5"(i?-2)(log(l+i?2))i-2" 

+ ^^^'^^f^5°(i?-^)(iog(i + i?^)r ^y 

Since H^l^ ^ is cubic at 0, we can pull out a factor of a and see that this 
belongs to 



which is admissible as we saw above. We consider -ffo(^'2fe) next. We write 



/(«) 



q;2 — sin a 
a"^ sm a 
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a^G^i^ lia) with G^^. ; E Q. There are two terms to consider. The first is 



and note that / is an analytic function of 63. We also write VF|j,;(a) 



X^sin^. ,, ,V V t'a'Gi,,{a) n .: 



a^ sm a 

3 



i=0 j=0 ^ ' \ J 

= /wEE (,.,).( l-,i (^ ■"«('+ «')) 



as desired. The other term is similar. Write 



sm a ~ a cos a 
1(a) = —. , 



and note that g is an analytic function of 63 . The remaining term in Hq is 



t^g{a)adaV2k- 



We expand V2k as before. If ada hits Gjj, ;, the fact that ada maps Q to Q' 
allows us to estimate the resulting term as above. Similarly ada{a^) — io?, so 
the contribution can be included in (1131) if ada hits a^ . It remains to consider 



a9„((ilog(l + i?^)))' = ^(ilog(l + i?^))'"\ 

which is even better than what we already had. This completes the analysis 
of HQ{v2k)- Finally, we consider the nonlinearity N2k{v2k)- This is dealt with 
exactly as in [T3], but we reproduce the calculations for completeness. Summing 
up the Vj we have 



U2k-i -uoe -^IS^iRlogR, Q). 

We will write /(a) for a generic even analytic function of a. Using lemma [5] 
in the appendix we get 
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t r^2 ^^2fc = /(a)^^(l - cos(2u2fe-i))w2fc 



sm a 



i?2 

iiW^ / 1 \2 



fe-i 

X 

fe-i 



E 77^7TW77^WiTI^^^'(^'(l°g^)""^ 2) 



(i'^A)2j(tA)2(fe+i-j) 



^ /55(i?(logi?)2^Q) 



(U) 



2' 



^ /5^(i?3(logi?)2'=-i,Q) 



(U)2 

For the quadratic term we have 

2Sm(2u2fe-l)/-, .0 XX 

t . „ (1 - COs(2w2fc)) 






2 sin^ a 



e i^(/s'(s->,a) + pjj/s'(siogR,a) 

, ^ /^■^(i?3(iogii')4fe-2^g) 



fe-1 

X 



E fi'^A)2jYa)2(fe-j) 1 2^ 



(i'^A)2j(iA)2(fe-j) V^^ (i''A)2'(iA)2('=+i-0 



^■^ (i-A)2'(a)2(fc+2-o 

2 



FinaUy, the cubic term can be bounded as 
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t T^ (2v2k - SU\{2v2k)) 

sm a 



{tX)^^'^IS^{RH\ogR)^^-\Q)\3 



R2 \z^^ ^iuxyjitxyik-j) 



fc-i 



c 



Z^ liu 



I S^W [log Rf^-'^.Q) 






£EE (,;;'V-^, "'(«(io8fi)->-,a) 



2 k-l 



^EE ft.,)..(a).(.-.) ^^^wlogi^)-^^go. 



This finishes the proof of Theorem [51 

3. The Perturbed Equation 

We wiU now start the construction of a complete solution from the approximate 
solution of the previous section. We let 

u(i, a) = U2k-iit, a) + e{t, a). 

This section is devoted to deriving the equation e has to satisfy for u to be 
a solution, and recasting this equation in a new coordinate system which allows 
easier analysis of the problem. In the a — t coordinates s has to satisfy 

cos(2g(Aa)) 

-£tt+£aa + cot aSa T^ £ ^ N2k-l[£) +62^-1, 

sm a 
or 

1 cos(2Q(Aq!)) ,^ , , ^^ , s ,„„x 

a a'' 

Here e2/c-i is the error at step 2/c — 1 from the previous section. N2k~i is 
defined in ([6]), but with U2k replaced by W2fc-i, and Hi is given by Q. We 
change variables to i? = X(t)a and t = v"^ "^ (so ^ = ■^), and let v{t,R) := 
£(i(r), A~^i?). Inserting this into 1^^ and rearranging we get 

dr + y ^5fl) ' + X (^- + ^^^^ 

= ^[iV2fc_i(£) + Hiie) + e2fc_i](i(T), A-ii?). 
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Our next goal is to transform dj^+^dn— '^"^^ ^i '' into a self adjoint operator 

on L^(R+, dR). To accomplish this we let £(r, R) :— R2v{t, R), and rewrite the 
equation above in terms of e as 



''"7->o \ ^ / ''" \2 I -^ o / '^^ 



-i^dr + -RdH) +^(-r+29ri-))e-Ce 



X-^R2{N2k-i{R'^s) + H,{R-2£) + e2k-i).. (30) 



where 



(|30p is the main equation we have to solve in this paper. C is self adjoint 
with respect to L^ (R^ , dR) (see the appendix for the definition of the domain of 
C) and its spectral properties are summarized in lemma S] of the appendix. We 
will adopt the notation of lemma |4] in the remainder of this paper. In particular 
recall the definition of the distorted Fourier transform associated with C 

■^ : / ^ /(O = lini / (t^{R,i)f{R)dR, 
with inverse 

-F-i : / ^ /(r) = hm T 0(i?, 0/(0^^)^6 

Here p{^)d^ is the spectral measure of C The idea is to expand e in terms of 
the generalized Fourier basis <j){R, ^) as 



/"OO 

e{T,R)^ x{T,R)cP{R,Op{Od^, 

Jo 



and derive a transport like equation for the coefficient x(r, 5) from pop . The 
operator ROr in ^U\i is not diagonal in the Fourier basis, but we can use the 
construction in |14j to deal with the error operator /C defined by 

Rd^J. = -2^d^u + K.U. (32) 

With K. defined as such, we have 





:F(dr 


+ 


>^] 


= (5r + y(-2eae + /c))j-. 


which yields 










F{dr + 


^us^y 


= {dr + 


^(-2eae + /c))V 








-(0,- 


^2^d^yT+2^lc[dr-- 








4(^ 


:2-^/C + 2[/C,C9d)-^- 



J^]^ 
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We can now apply our Fourier transform to (j30p to get (with x = J-e) 



■'''^^■'-\e,{^)U 



(33) 

Our strategy will be to solve this equation for the Fourier coefficient x and 
use the inverse Fourier transform T^^ to retrieve a solution of (1511)) . Note that 
we are interested in solutions of (1331) which decay as r — )■ oo . This means that we 
have to solve the equation backwards in time with zero Cauchy data at t = cxd. 
The next section is devoted to describing the mapping properties of the error 
operator /C and the fundamental solution of the "transport like" equation 






x(T,i) = h[T,i). (34) 



We call this a "transport like" equation because the characteristic curves of 
the operator 9r — 2^^9^ are (r, A~^(t)^) (see section 8 of [2] for more details). 



4. The Fundamental Solution and Estimates in H" 

In this section we define certain Sobolev spaces associated with the operator £. 
We will also provide the mapping properties of the fundamental solution opera- 
tor of p4p and the operator K. relative to these spaces. Finally we demonstrate 
how to control the norm of the right hand side of p3p . In the next section we will 
use a contraction mapping argument on these Sobolev spaces to find a solution 
£ to (jSni) with the right decay properties. The estimates from this section are 
the necessary ingredients for being able to close that fixed point argument. 

We begin with the norms. On the frequency side we define the weighted L^ 
norms 

ii/iL-:=(y^ l/(OHo'v(e)rf?)^ (35) 

For functions of the spacial variable R we define 



The relationship between this norm and the usual Sobolev norm on M^ is 
explained in lemma [B] in the appendix. Finally to control the decay in time we 
introduce the L^'°°Lp''' spaces with norm 
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||/|li=o.iVi2,. :=supT^||/||^2,= . 

L^'°°Hp is defined similarly. 

We now move to the estimates. A priori we only have 

/C:Co"((0,(X)))^C°°((0,(X))). 

However, with the norms just defined, according to lemma[S]in the appendix, 
/C has the following mapping properties 

1C:LY^lY^^, (36) 

[K.,£,d{^:LY^LY. (37) 

Next let H be the backward fundamental solution of the operator 

5.-2^eae)'+e, (38) 

and denote its operator kernel by H[t, a). By this we mean that suppressing 
the ^ variable, 

is a solution of ([M|) . Then according to lemma [7] in the appendix, given s 
there is constant C — C{s) such that for large enough N 



'ioc,JV-2i 



A 



,W-1^2 



Motivated by this we want to bound the L°°-'^ L?'^ norm of the right hand 

side of dMl) by the L°°'^-^lI''^^ norm of x and the L^^^^-^L^.^ norm of {dr - 
2^S,d^)x. The terms involving K. or its commutator can be bounded using p6p 
and (|37p . The term coming from e2k-i will be dealt with in the next section. 
Note that this term does not depend on x so we only need to make sure that 
it belongs to our iteration space (to be specified in the next section). This will 
be accomplished by taking k large as in [U:. For the N2k~i term we will use 
lemma [5] from the appendix. Therefore the term with Hi is the only one which 
is essentially new. The estimate we need is 



Recall that 



||A-^i?5ij^(i?-2£)|| „ < ||e|| 1. (40) 



sin a — a cos a sin a — a^ 

Hi{v) = : daV H -^-^ 1 

a sm a a^ sm a 
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We write 



so 



sm a — a cos a 

z(a) = 5-^ , 



V a^ sin a ' 

Since A^^ = 07^"^^'^ and according to lemma|ni the first term above can be 
bounded consistently with (|40)) . For the second term note that as long b& v <\ 
we have A^'^A = ^-"^+1+2!^ = cT~'^t^~^ < r~^, which is he desired time gain in 
(Uni)- Note moreover, that since we are only concerned with finding a solution 
inside the light cone, we may replace r by q{^)^ where g is a smooth function 
supported in i? < 1. Therefore, again revoking lemma O for such range of v it 
suffices to bound 



R2q{ 



R\ 



RdRiR~-^e)\\Hi<\\el^.+ i. 



A " H_ 

But noting that in two dimensions, ROr — xdx + ydy, lemma |6] allows us to 
reduce this to 



A' 

kilo Vn^li^o -J-T-nd V\iar>omciQ l-\dr>ono<ii _^ ^11.1 



where f — R ^e. This holds true because because ^ and ^ , together with 



all of their derivative, are bounded (see e.g. [iU]). 



5. The Blow Up Construction: Proof of Theorem [T] 

In this section we combine the results from the previous sections to find e in a 
way that leads to blow up for u. Fix ly S (^ 1] and let U2k-i and e2k~i be as in 
theorem [21 The index k is chosen sufficiently large depending v and in particular 
on N (as will be specified below). A priori U2fc-i and e2k-i are defined only on 
the cone {a < t}, but we extend them to functions on the double cone {a < 2t}, 
with the only requirement that extensions are of the same smoothness class and 
all meaningful derivatives agree on the boundary of the cone. With these choices 
of ly, U2k-i and e2k-i, and with s G (|, ^) we will find a solution e of (|30p (solved 
backwards in r) which satisfies 

||£(r)||^.,. <r^-^, ||(a. + ^i?a^)e(r)||^^ <r^-^. (41) 

The large exponent TV is determined by lemmas [7] and [51 In what follows we 
identify the Sobolev spaces on the domain sphere (more precisely on a neighbor- 
hood of the north pole) with the normal Sobolev spaces on K^ via the {a, 9) (or 
{R, 9)) coordinates. 
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Near the boundary of the cone {a = t}, the error 62^-1 has a singularity of 
the type (1 — aY~^ log™(l — a), which means that locahy near the boundary 
e2k-i G H^ as long as /3 < z/. Note that under the transformation T of lemmajS] 
this corresponds to Hp . On the other hand according to theorem [5] 

Near R — this means that we have an expansion of the form (with T as in 
lemma [5]) 

T{Rh2k~i) - e'^i?(co(T) + ci(T)i?2 + C2(r)i?4 + ...), 

which is smooth around R — 0. Finally we consider the size of the error for 
large R 

•i?(log(2 + i?))2^-i 
e2fc-i = U' 



t-2fe+2i2(iA)2 

It follows that the L"^ norm of T(A~2i^2 62^-1) is bounded by CT'^t^''~'^ (note 
that we are taking the L^ norm only on the light cone which corresponds to 
R^t). Taking k so large that 

supr^+2i2/c-2<^^ 

r>l 

and noting that due to the expansion of 62^-1 for large R taking R derivatives 
reduces the size for large R, the arguments above show that 

\\\-^R--e2k-i\\L^.^H'^<l, (42) 

as long as s < -^ . Now using the distorted Fourier transform F we recast (j30l) 
as ((331) with x = Fe. According to lemmas [3 H [3 H and (gO]) and (|42|) we can 
solve p3p using a contraction mapping argument with respect to the norm 



1^".= 



Ar, 



(dr - 2^^d^ 



l^-^-^lI-" 



Using the inverse of the distorted Fourier transform and lemma [5] we can 
derive a solution e of (l30ll which satisfies (l4Tt. 



To go from the longitudinal variable m to a full co-rotational wave map in 
terms of the ambient coordinates of R'^ 3 5'2, observe that these coordinates are 
given by (/) o T(u), where (j) -."^^ ^ S'^ C'M? \s given by 



(j){e v) ~ (cos ti, sin w cos 6*, sin w sin 0). 

It is then verified that (j)oT{u) G i/2«+i^ interpreted componentwise. We have 
thus constructed a wave map on the closure of the cone {a < t, < t < to}, 
which is of class H^'^'^~ . To extend this to a wave map on all of 5'2 x R, extend the 
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initial data dtu(to, •), u{tQ, •) at time t = to to all of S^ in the same smoothness 
and equivariance class, and in such a way that a large neighborhood of the 
south pole (a — tt) is mapped to the south pole on the target. Let ii be the 
corresponding wave map. We claim that u does not develop singularities in 
(0,io] X 5^. Indeed due to the equivariance assumption the first singularity has 
to develop at one of the poles. At the north pole a = this is precluded by our 
construction because finite speed of propagation implies that u agrees with u on 
the light cone {a < t, < i < io}- At the south pole this is precluded by the 
small energy result in [35] which implies that singularity development results in 
energy concentration, which in turn is ruled out by the constancy of the initial 
data near this pole and finite speed of propagation. This completes the proof of 
theorem [1] 



6. Appendix 

In this section we provide the statements, without proofs, of some of the results 
from .14. which are used in the current work. 

Lemma 1. Let fc > 1, and let L be defined as in jj^p . Then the solution v to 
the equation 

Lv^fe S\R-^ {log Rf''-^), v{0) = v'{0) = 

has regularity 

ve S^{R{\ogRf''-^). 

Proof. Lemma 3.7, p. 560 in [14]. 

Lemma 2. Let 

1 



ve--^LS^RlogR,Q). 



The' 



n 



sint>e -^LS^{R\ogR,Q), cosv e LS°{1,Q). 

{txy 

Proof. Lemma 3.8, pp. 561-562 in [14| . 

Lemma 3. Let f be an analytic function in [0, 1) with an odd expansion at 0, 
and let (3 > -^ be a fixed real number. Then there is a unique solution w to the 
equation 

Lpw - /, w{Q) = A, daw[Q) = B, 
with 

Lfi = {l- afdl + (fl-i + 2ali - 2a)da + (-/9^ + /? - a^), 
such that 
(i) if A = B = 0, then w is analytic in [0, 1) with an odd cubic expansion at 0. 
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(it) 

w{a) =ci0i(a) + C202(a) 

-(/3 + l/2)-Vi(a) / ^2{a')qi{a)I{a)da' 



(/3 + l/2)-V2(a) / 0i(a')gi(a')/(«')rfa'- 



Here ao < 1 is some number close to one, and ci and ci are constants determined 
by the initial values, qi is given by 

q,{a) = {l-a)-f'-^l + {\-a)Ua)] 

with 92 analytic near a — \. Moreover we can find constant fj,i and jli such 
that 



</.i(a) = 1 + ^ Mi(l - a)\ Ma) = (1 - «)^+^ [l + E^~'(l " ")' 
1=1 1=1 

'i'f P ~ ^ ^ ^^ ■ If (3 ^ ^ ^ ^^ then we need to modify <j)i to 

oo 

01 (a) = l + ^fii{l -a)' +C(/)2(a)log(l-a), 
1=1 
for some constant c depending on (3. 

Proof. All statements come from lemma 3.9 in [14]. Part (i) corresponds to part 
(i) there. Part (ii) comes from equation (3.31) in the proof of lemma 3.9 in [H] . 
The expressions for (f)j are given by (3.27) and (3.30), and for qi at the bottom 
of page 567, in 14 . 

Lemma 4. Let C be the self adjoint operator —d^A-^^—ji^^Tyi on L^{{0,oo),dR) 
with domain 

Dom(£) = {/ e L^({0,^))\fj' e ACioc{(0.^)),Cof E L^((0,^))}, 
where Cq := —d\ + j^. Then: 

(i) The spectrum of C is purely absolutely continuous and equals spec(C) = 

[0,oo). 

(ii) For each z G C there exists a fundamental system (j){R,z), 9{R,z) for 
C — z which is analytic in z for each R > and has the asymptotic behavior 

(j){R, z) r^ Ri , e{R,z) r^ -R~^ asR-^O. 

In particular, their Wronksian is W{9{., z),(j)(., z)) — 1 for all z G C. By 
convention 9{R, z) and (j){R, z) are real valued for z G M. is the Weyl- 
Titchmarsh solution of C — z at R — 0, and admits an absolutely convergent 
asymptotic expansion 
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(j){R, z) = MR) + R'^ Y^^R^^y^iiR^) 



whe 



MR) 



Ri 



l + i?2' 



and for j > 1 the functions (j)j are holomorphic in U — {Re m > — 2} 'J'^'^ 
satisfy the bounds 



\<P,{u)\<j--^\og{l + \u\). 



In particular 4>j{0) = and \(j)'JO)\ < (^_i)[ for all j 



> 1. 



(Hi) For each z G C, Im z > 0, let ip^{R,z) denote the Weyl-Titchmarsh 
solution of C ^ z at R = 00 normalized so that 



^p'^(R,z)^z ie'""^^ 'AsR^oo, Im z^ > 0. 

If^ > 0, then the limit ip^ (R, ^+iO) exists point-wise for all R > and we de- 
note it byip+{R,^). Moreover, define ip^ {., ^) := V'+(-,C)- Themp+{R,^), ■0"(i?,C) 
form a fundamental system of C — ^ with asymptotic behavior tp^{R,(_) ^ 



C 



-i ±J«5fl 



as i? —> 00. More precisely, for and ^ > we can write 



where a admits the asymptotic series approximation 






with zero order symbols ipj' that are analytic at infinity, 



svLp\{RdR)''ij+{R)\<^ 



R>0 



in the sense that for all large integers jo, and all indices a, {3 we have 



sup 

K>0 



{RdRriqd.f [aiq, i?) - ^ g~-'V'+(i?) 



< Ca, 13,301 



-jo-1 



for all q > 1. 
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(iv) The spectral measure p{S.)d^ of C is absolutely continuous and its density 
is given by 

PiO = -Im m{£, + iO)x[^>o] 
with the "generalized Weyl-Titchmarsh" function 

W{9i.,zUH.,z)) 
"^^'^ VF(^+(.,z),0(.,z))' '™"-"- 



Moreover p satisfies □ 



and we can write 



p{0^{f°^'^' \tl 



p{0 = -HO\-^ (43) 

TT 

where a satisfies the symbol type bounds 

\mfa{o\<ck\am ve>o. 

(v) The distorted Fourier transform defined as 

■^:/^/(0= lim / <P{R,Of{R)dR 
is a unitary operator from L^(R"*') to L'^{M.'^,p), and its inverse is given by 

T-^:f^f{r)= lim T 0(i?, O/COKCK- 

Here lim refers to the L^(M+,/3), respectively the L^(R+), limit. 

Proof. Unless otherwise stated, "lemma x," "proposition x" or "theorem x" refer 
to the corresponding lemma, proposition or theorem in [Tl]. Part (i) is lemma 
5.2. Part (ii) is part (a) of theorem 5.3. For the asymptotic behavior of see 
lemma 5.4. Part (iii) comes from part (b) of theorem 5.3 and proposition 5.6. 
(iv) comes from (c) of theorem 5.3 and proposition 5.7. (v) is (d) of theorem 5.3. 
Also see f5^, and in particular section 3 and example 3.10 there, for a general 
analysis of Schrodinger operators with a singular potential similar to C 



Lemma 5. The operator K from i32\) maps 
In addition we have the commutator bound 
Both statements hold for all s € R. 



means that for some constant C one has C ^a < b < Ca. 
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Proof. This is the content of proposition 6.2 in [Tl]. See also theorem 6.1 there. 
Lemma 6. Define the map 

u{R) ^ Tu{R, 9) = e''^R-^u{R), 

where the right hand side is interpreted as a function in M? expressed in polar 
coordinates {R,9). Then T is an isometry 

T ■.L^{R+) -^ L^{R^), 
and for any s > we have 



||u||^./.(R+)x||ru| 



H-( 



in the sense that if one side is finite then the other is finite and they have 
comparable sizes. 

Proof. See lemma 10.1 in [T3]. 

Lemma 7. Let H be the fundamental solution of the operator i38\) . Given s > 
let N be large enough. Then there is a constant C which depends on s but not 
on N so that 



i^^iU «.,.^.=+i + II (^- - 2xea«)mL ._, r... ^ cN'^m^^^.^. 



'L°°."-2i^ 



L°°.«-iL?' = 



Proof. See proposition 7.1 and corollary 7.2 in |14) . 

Lemma 8. Assume that N is large enough and that ^ < s < f + f . Then the 
map 

X ^ X-^T{R^ N2k-i{R^^ T^^x)) 

is locally Lipschitz from L°°''^-^lI''^^ to L°^^^ Ll^' . 
Proof. Proposition 7.3 in |14) . 
Lemma 9. Let q e S{1, Q) and \s\ < ^ + |. Then 

Ik/ll//^ < ll/IU- 

Proof. This is lemma 9.1 in [14]. Note that the algebras in our paper are slightly 
different from those in |14j . In particular there is an extra dependency on the 
variables 62, • • • ,65. This extra dependency is irrelevant in the proof of the lemma 
and can be ignored just as the 61 variable is ignored in the proof of the lemma 
in tl4i. 
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